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ABSTRACT 

This paper summarises an investigation of the effects of low amplitude noise and pe- 
riodic driving on phase space transport in three-dimensional Hamiltonian systems, a 
problem directly applicable to systems like galaxies, where such perturbations reflect 
internal irregularities and/or a surrounding environment. A new diagnostic tool is ex- 
ploited to quantify the extent to which, over long times, different segments of the same 
chaotic orbit evolved in the absence of such perturbations can exhibit very different 
amount of chaos. First passage times experiments are used to study how small per- 
turbations of an individual orbit can dramatically accelerate phase space transport, 
allowing 'sticky' chaotic orbits trapped near regular islands to become unstuck on 
surprisingly short time scales. The effects of small perturbations are also studied in 
the context of orbit ensembles with the aim of understanding how such irregularities 
can increase the efficacy of chaotic mixing. For both noise and periodic driving, the 
effect of the perturbation scales roughly logarithmically in amplitude. For white noise, 
the details are unimportant: additive and multiplicative noise tend to have similar 
effects and the presence or absence of a friction related to the noise by a Fluctuation- 
Dissipation Theorem is largely irrelevant. Allowing for coloured noise can significantly 
decrease the efficacy of the perturbation, but only when the autocorrelation time, 
which vanishes for white noise, becomes so large that there is little power at frequen- 
cies comparable to the natural frequencies of the unperturbed orbit. This suggests 
strongly that noise-induced extrinsic diffusion, like modulational diffusion associated 
with periodic driving, is a resonance phenomenon. Potential implications for galaxies 
are discussed. 

Key words: chaos - galaxies: formation - galaxies: kinematics and dynamics 



1 MOTIVATION 

The objective of the work described here was to explore 
phase space transport in complex Hamiltonian systems that 
admit both regular and chaotic orbits and, especially, to 
understand how and why low amplitude perturbations, ide- 
alised as noise and/or periodic driving, can dramatically ac- 
celerate the rate at which a single chaotic orbit moves from 
one part of phase space to another. 

This work has immediate implications for a variety of 
problems related to galactic astronomy. In a first approxi- 
mation, an elliptical galaxy can perhaps be characterised by 
a smooth, time-independent bulk potential. However, this is 
surely not the whole story. One must also allow for discrete 
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substructures, e.g., individual stars, and the presence of a 
surrounding environment, both of which can have apprecia- 
ble effects under appropriate circumstances. 

It has been long accepted that discreteness effects, i.e., 
gravitational Rutherford scattering between nearby stars, 
can be modeled as white noise and friction by a Fokker- 
Planck equation. Successive encounters are idealised as a se- 
quence of instantaneous kicks associated with random forces 
that are delta-correlated in time, so that the forces at two 
different instants are statistically independent. This noise is 
then augmented by a dynamical friction which represents the 
systematic drag associated with a star moving through the 
ambient medium. In the context of Chandrasekhar's (1943a) 
binary encounter approximation, the amplitudes of the fric- 
tion and noise are connected by a Fluctuation-Dissipation 
Theorem (Chandrasekhar 1943b), a generic result to be ex- 
pected in any statistical treatment of an isolated system (cf. 
van Kampen 1981). 

Regular motions associated with satellite galaxies and 
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other companion objects can trigger a near-periodic driving 
that will also effect stars in the parent galaxy. In terms of 
the parent mass M, the companion mass m, the parent size 
r, and the characteristic separation d between the galaxies, 
a typical amplitude and frequency are easily estimated. De- 
noting, respectively, by v and u the typical speed of a star 
in the parent and the relative speed of the parent and com- 
panion objects, the perturbing force will have an amplitude 
of order (m/M)(r /d) 2 times as large as the force resulting 
from the bulk potential of the parent galaxy, and the char- 
acteristic frequency lo will scale as lo ~ (r /d)(u/v)t~j D 1 , with 
io a characteristic dynamical time for the parent galaxy. 

A galaxy situated in a high density environment, e.g., 
a rich cluster, will feel a superposition of forces from other 
nearby galaxies which, generically, is far from periodic (al- 
though there could be a near-periodic component if, e.g., the 
galaxy has bound satellites). It seems reasonable to suppose 
that these forces are 'random,' but it is obviously not rea- 
sonable to pretend that they correspond to instantaneous 
events. Idealising them as white noise is clearly inappropri- 
ate. However, following standard techniques from statistical 
physics, what does seem appropriate is to model them as 
coloured noise (cf. Honerkamp 1994), allowing for successive 
impulses which, albeit random, have finite duration. 

Assuming that the noise is a Gaussian random process 
with zero mean, it is characterised completely by its second 
moment, which is assumed to satisfy 

(F a (ti)Fb(t2)) = K(ti-t 2 )8ab (a,b = x,y,z), (1) 

with K the autocorrelation function. Demanding that K(t) 
be proportional to a Dirac delta Sd(t) yields white noise, 
with a flat power spectrum and vanishing autocorrelation 
time. Allowing K to be nonvanishing for a range in r 
yields coloured noise, with a band-limited power spectrum 
and a finite autocorrelation time. If, consistent with Chan- 
drasekhar's (1941) nearest neighbor approximation, one as- 
sumes that irregularities associated with the presence of 
other galaxies reflect primarily the effects of a few nearby 
neighbors, a characteristic amplitude and autocorrelation 
time are again easily estimated. 

Given that the typical distance between galaxies is often 
less than ten times the size of a typical galaxy, an external 
environment can easily induce perturbations with amplitude 
at the one percent level or above. The power associated with 
these perturbations should peak at frequencies somewhat 
smaller than . 

But why might such weak perturbations matter? For ex- 
ample, conventional wisdom holds that discreteness effects 
should be completely irrelevant on time scales short com- 
pared with the binary relaxation time tn, a time which, for 
a galaxy like the Milky Way, is orders of magnitude longer 
than the age of the Universe tn ■ The answer here lies in the 
fact that galaxies are more complicated dynamical systems 
than has been appreciated until recently. High resolution 
photometry has provided clear evidence that many galaxies 
are genuinely triaxial, i.e., neither spherical nor axisymmet- 
ric, and that most galaxies probably have a high density 
central cusp. Indeed, these features are so well established 
that they have been proposed as the basis of a new classifi- 
cation scheme for ellipticals (Kormendy and Bender 1996). 
However, dynamical considerations involving the manifes- 
tations of resonance overlap provide compelling reasons to 



believe that the combination of cusps and triaxiality leads 
generically to a bulk potential corresponding to a complex 
phase space that admits significant measures of both regular 
and chaotic orbits (cf. Merritt 1996, 1999) (although there 
are simple examples of cuspy triaxial potentials which are 
completely integrable [cf. Sridhar and Touma 1996,1997]). 

A chaotic bulk potential can have profound implications 
for the behaviour of stars moving in a galaxy. The sensitive 
dependence on initial conditions characteristic of individual 
chaotic orbits (cf. Lichtenberg and Lieberman 1992) implies 
that initially localised ensembles of chaotic orbits tend to di- 
verge exponentially, which leads to a phase mixing far more 
efficient than what obtains for ensembles of regular orbits 
(Kandrup & Mahon 1994a, Mahon et al 1995, Merritt & 
Valluri 1996, Kandrup 1998). This chaotic mixing has po- 
tentially significant implications for the efficacy with which 
such irregularities as metallicity gradients can disperse in a 
near-equilibrium galaxy. To the extent that galaxies out of 
equilibrium are dominated by chaotic orbits, chaotic mixing 
could also provide a compelling explanation of the obvious 
efficiency of violent relaxation (Lynden-Bell 1967), as ob- 
served, e.g., in numerical simulations. 

However, the complex character of the phase space as- 
sociated with a potential admitting both regular and chaotic 
orbits significantly limits the efficiency of chaotic mixing in 
a time-independent potential. Although the chaotic phase 
space on a constant energy hypersurface is usually con- 
nected in the sense that, in principle, a single orbit can 
and will eventually access all of it, partial obstructions like 
cantori (cf. Percival 1983) in two-dimensional systems and 
Arnold webs (Arnold 1964) in three-dimensional systems can 
severely restrict phase space transport. Although an orbit 
eventually passes from one phase space region to another, 
it can only do so by traversing bottlenecks which so impede 
progress that the time scale for the transit can be thou- 
sands of dynamical times, this corresponding to a time much 
longer than tn- Indeed, for three-dimensional systems one 
anticipates that, in many cases, the time to diffuse along an 
Arnold web is exponentially long (Nekhoroshev 1977). One 
manifestation of this fact is the so-called 'stickiness' phe- 
nomenon noted by Contopoulos (1971), whereby a chaotic 
orbit can become stuck near a regular island and behave in 
a near-regular fashion for a time t 3> tn- 

The crucial point is that even very low amplitude per- 
turbations can dramatically accelerate phase space trans- 
port through these bottlenecks, decreasing the time required 
for an orbit to transit from one phase space region to another 
and, consequently, the time required for an initially localised 
ensemble to probe all the accessible phase space. This was 
first recognised in the context of the so-called Fermi accel- 
eration map (Fermi 1949), a simple symplectic map used 
to model cosmic ray acceleration, where the introduction of 
white noise was shown to greatly accelerate the rate of phase 
space transport (Lieberman & Lichtenberg 1972). More re- 
cently, Habib, Kandrup, & Mahon (1996, 1997) revisited 
the role of white noise in phase space transport, which as- 
sumes practical importance to accelerator dynamicists con- 
cerned with high density charged-beam experiments, where 
discreteness effects can result in the degradation of an ini- 
tially focused beam (Habib & Ryne 1995). Accelerator dy- 
namicists have also recognised that periodic driving can act 
comparably as a source of accelerated phase space trans- 
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port. If orbits are pulsed with a frequency comparable to 
their natural frequency, the resulting resonant coupling can 
dramatically alter the rate at which the orbits move through 
phase space (Lichtenberg & Lieberman 1992). 

Low amplitude perturbations may be especially impor- 
tant in near-Hamiltonian systems like galaxies where the 
bulk potential is generated self-consistently by the stars 
themselves, rather than being imposed externally. Dating 
back at least to Binney (1978), galactic dynamicists have 
assumed that nontrivial structures like triaxiality in an el- 
liptical or the bar of a spiral cannot be supported completely 
by chaotic orbits. Rather, such interesting structures would 
seem to require the presence of various regular orbit families 
which coexist with the chaotic orbits in the complex poten- 
tial and which, owing to their shape, can serve as a skeleton. 
However, several groups interested in modeling galaxies have 
found that, because of resonance overlap, there may exist al- 
most no regular orbits in crucial phase space regions (e.g., 
near corototation), which might suggest that self-consistent 
models do not exist. One proposed solution (cf. Athanas- 
soula et al 1983, Wozniak 1994) has been to build the skele- 
ton with 'sticky' orbits which, albeit chaotic, can behave 
as if they are nearly regular for times long compared with 
the age of the Universe. Indeed, in the context of numeri- 
cal modeling based on Schwarzschild's (1979) method, both 
Merritt & Fridman (1996) and Siopis (1998) have concluded 
that near-equilibria appropriate for triaxial Dehnen poten- 
tials cannot be constructed without including a substantial 
number of nearly regular chaotic segments. 

This seems completely reasonable, but one needs to be 
sure that the resulting models are sufficiently stable towards 
small perturbations reflecting short scale irregularities in 
the potential or the effects of nearby galaxies. The problem 
then is that 'sticky' orbits tend to become unstuck much 
more quickly in the presence of low amplitude irregulari- 
ties than they would in the absence of such irregularities. In 
particular, numerical experiments (Habib, Kandrup, & Ma- 
hon 1997) have suggested that extremely weak white noise 
perturbations corresponding to a relaxation time as long as 
t R ~ 10 B - W 9 t D can have major effects within ~ lOOto, a 
time shorter than the Hubble time, even if an unperturbed 
flow remains unchanged for ~ 500to or more. Consistent 
with Merritt (1996), one might therefore anticipate that, in 
response to low amplitude irregularities, galaxies will tend 
systematically to evolve from strongly triaxial configurations 
towards more nearly axisymmetric configurations. 

Given the recognition that low amplitude perturbations 
may play an important role in the structure and evolution 
of galaxies, there are three obvious questions to address: 

1. What is the physics that triggers accelerated phase space 
transport? If nothing else, understanding why this phe- 
nomenon arises should help one develop an intuition as to 
when it might prove important. 

2. How does the effect scale with amplitude? Answering this 
should enable one to estimate dimensionally when it is that 
these effects must be considered. 

3. To what extent do the details matter? Will one get wildly 
different behavior for additive (i.e., state-independent) and 
multiplicative (i.e., state-dependent) noises, or are such 
state-dependent effects largely irrelevant? Under what cir- 
cumstances does allowing for colour, i.e., endowing the ran- 
dom forces with a finite autocorrelation time, alter the effi- 



cacy of the noise? To the extent that the details are unim- 
portant one will have the luxury of being able to ignore many 
complications which are nearly inaccessible observationally. 

The computations described below suggest strongly 
that, like periodic driving, noise-induced phase space trans- 
port is a resonance phenomenon which involves a coupling 
between the frequencies at which the noise has substantial 
power and the natural frequencies of the unperturbed or- 
bit. It follows that allowing for a finite autocorrelation time 
t c has only a minimal effect provided that t c remains short 
compared with tn, the natural time scale for the unper- 
turbed orbits. However, as t c increases the efficacy of the 
force decreases and, for t c 2> to, the effects of the noise be- 
come relatively minimal. Overall, the effects scale logarith- 
mically in amplitude. The dependence on t c is more com- 
plex, but is again better represented as logarithmic than 
any simple power law. Other aspects of the perturbation 
seem largely immaterial. Additive and multiplicative noises 
tend to have virtually identical effects, and the presence or 
absence of dynamical friction does not seem to matter appre- 
ciably. In this sense, this problem of diffusion through bot- 
tlenecks, so-called entropy barriers (cf. Machta & Zwanzig 
1983), is very different from energy barrier problems, where 
the form of the noise tends to matter a very great deal (cf. 
Lindenberg & Seshadri 1981, Alexander & Habib 1994). 

In the past, a number of individuals, including the first 
author, have suggested that noise is important as a source of 
accelerated phase space diffusion primarily because it allows 
one to violate the Hamiltonian constraints associated with 
Liouville's Theorem, which makes hunting through phase 
space easier. This does not seem to be correct. The physics 
seems to be essentially the same as for periodic driving, so 
it is clearly not the non-Hamiltonian character of the noise 
per se that is responsible for what is seen. 

Section II discusses some basic issues related to phase 
space transport in complex time-independent potentials in 
the absence of all perturbations. Section III studies the ef- 
fects of low amplitude perturbations on individual orbits by 
performing first passage time experiments. What this entails 
is selecting an orbit which is originally trapped near a reg- 
ular island and determining how the characteristic escape 
time depends on the amplitude and form of the perturba- 
tion. Section IV focuses on how chaotic mixing is altered 
through the introduction of noise and periodic driving. Sec- 
tion V concludes by discussing potential implications for real 
galaxies. 



2 PHASE SPACE TRANSPORT IN CHAOTIC 
HAMILTONIAN SYSTEMS 

The computations described here were performed for three- 
degree-of-freedom Hamiltonian system of the form 

H= l(pl+pl+pl) +V(x,y,z), (2) 

with V given as a generalisation of the two-degree-of- 
freedom dihedral potential (Armbruster et al 1989), with 
two free parameters, a and b: 

V(x,y,z) = -(x 2 +y 2 + z 2 ) + ^(x 2 + y 2 + z 2 ) 2 
-\{x 2 V 2 + ay 2 z 2 + bz 2 x 2 ), (3) 
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The two-dimensional limit of (3) appropriate for motion 
with 2 = served as a prototypical example in several ear- 
lier papers (Mahon et al 1995, Habib, Kandrup, and Mahon 
1997) which discuss its physical properties extensively. The 
fully three-dimensional version was explored in the context 
of chaotic mixing in Kandrup (1998). Note that, for the en- 
ergies used here, 1 < E < 20 or so, a dynamical time tr> 
corresponds to t ~ 2 — 5, and that most of the power in 
individual orbits is for uo ~ 1 — 5. 

It is well understood that, because of cantori, chaotic 
orbit segments in two-degree-of-freedom Hamiltonian sys- 
tems often decompose naturally into two (or more) distinct 
classes, namely (1) unconfined segments which look wildly 
chaotic and tend to avoid regions near regular islands, and 
(2) confined, or sticky, segments which are trapped near reg- 
ular islands and are nearly indistinguishable visually from 
regular orbits. Whether one should expect comparable dis- 
tinctions in three-degree-of-freedom systems is not com- 
pletely obvious. General topological arguments imply that 
a generic chaotic phase space will involve separate regions 
connected by an Arnold web (Arnold 1964), but there is no 
guarantee that different regions will be particularly near- 
regular or wildly chaotic. Nevertheless, it is clear from vi- 
sual inspection that a single chaotic orbit often does divide 
naturally into nearly regular and wildly chaotic pieces. The 
obvious question, therefore, is how to characterise such dis- 
tinctions in a robust and quantifiable fashion. 

One way is to compute estimates of the largest short 
time Lyapunov exponent (cf. Grassberger et al 1988) for 
different segments of the same chaotic orbit, and determine 
whether iV[x], the distribution of short time Lyapunov expo- 
nents, exhibits behavior suggestive of multiple populations. 
Such computations were effected here in the usual way by 
introducing a small initial perturbation, evolving the un- 
perturbed and perturbed initial conditions, and periodically 
renormalising the perturbed orbit to ensure that the pertur- 
bation remains small (cf. Lichtenberg & Lieberman 1992). 
Explicitly, this involved computing an approximation to 



1, \SZ(t) 
X(t) = lim - In ' ; ; 
Ay ' *z(o)^o t \SZ(0)\ 



(4) 



with \SZ\ 2 = \Sr\ 2 + \5p\ 2 . Note that, given x{ti) and x(^) 
evaluated in this way, the average exponential instability 
for the interval ti < t < satisfies (cf. Kandrup & Mahon 
1994b) 
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ta-ti 
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This tact works because chaotic segments which look 
nearly regular and are situated relatively close to regular 
islands tend systematically to have smaller short time Lya- 
punov exponents than more wildly chaotic segments. (This 
fact, which seems intuitive physically, can be quantified us- 
ing the notion of 'orbital complexity' [Kandrup et al 1997], 
which characterises the extent to which an orbit segment 
has considerable power in a large number of Fourier modes. 
Nearly regular chaotic segments tend to have low complex- 
ity, i.e., power concentrated in a small number of modes, but 
they also tend to have smaller short time Lyapunov expo- 
nents. If short time Lyapunov exponents % an d complexities 
n are computed for different segments of the same chaotic 



orbit, one finds typically that, for reasonable sampling times, 
the rank correlation 1Z{x,n) ~ 0.75 — 0.95.) 

For a large number of different potentials and energies, 
data were generated by selecting 16 different initial condi- 
tions in the same connected phase space region and integrat- 
ing each for a total time t = 524288 = 2 19 , with x(f) an d 
the phase space coordinates recorded at intervals At = 1. 
The resulting orbits were partitioned into segments of length 
t — 2 fc , with k a positive integer, and short time exponents 
determined for each segment using eq. (5). These short time 
exponents were then binned to extract distributions of short 
time Lyapunov exponents, JV[x]j an d the forms of these dis- 
tributions analysed as functions of k or \ogt. 

Suppose that, for some potential and energy, there is 
only one 'type' of chaotic orbit, and that the time scale on 
which the local instability of the orbit changes significantly, 
i.e., the autocorrelation time r for the 'local stretching num- 
bers,' is short compared with the times over which the orbit 
is being probed. Let St denote some basic interval and let 
J\f(5t) denote the distribution of short time Lyapunov expo- 
nents for an interval of this length. Assuming only that the 
moments of M[x(St)] exist, the Central Limits Theorem (cf. 
Chandrasekhar 1943b) then makes two specific predictions 
about longer intervals mSt for m 3> 1: 

(1) The longer time distribution N[x(t)] for t = mSt will be 
approximately Gaussian. 

(2) The relative width of this Gaussian will decrease as 
m _1//2 , so that the dispersion a x oc t -1 / 2 . 

Deviations from a Gaussian and/or a dispersion that de- 
creases more slowly are prima facia evidence that, over the 
time scale of in question, the orbit segments decompose into 
more than one distinct population. 

So what is actually observed? For short times, t < 
16 — 32 or so, one generally sees a singly peaked distribution 
iV[x]. However, this does not necessarily imply that there is 
only one population. For short t the distribution N[x] is so 
broad that one could easily miss a good deal of structure; 
and it is hard to exclude the possibility of several popula- 
tions with distinct peaks which, when convolved, still yield 
a unimodal distribution. As illustrated in the top panels of 
FIGS. 1 and 2, for small t the distribution is typically very 
smooth and resembles a Gaussian. However, there usually 
are statistically significant differences from a normal distri- 
bution. (Panels (a) in FIG. 1 and 2 were each generated by 
binning 2 23 » 8.4 x 10 6 segments, so that even small irreg- 
ularities are significant!) In particular, N[x] typically has a 
pronounced skew, as might be expected generically if the to- 
tal distribution is comprised of several distinct populations. 

As t increases, one sometimes sees evidence for multiple 
populations, but not always. The evidence for multiple pop- 
ulations arises invariably for those potentials and energies 
where visual inspection indicates the possibility of trapping 
near a regular region for reasonably long periods of time. 
When such trapping is not observed, for t > 32 or so N[x] 
typically corresponds very closely to a true Gaussian distri- 
bution, with no statistically significant skew. The smallest 
values of x ar e usually significantly larger than zero. Alter- 
natively, when trapping occurs with reasonable frequency 
N[x] tends instead to reflect the sum of a Gaussian distri- 
bution centered at a comparatively high value of x plus one 
(or more) additional lower-x populations which can be man- 
ifested as contributing a secondary peak to N[x] and/or an 
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extended low-x tail extending down to very small values of 
X- Examples of potentials and energies which exhibit only 
one and more than one populations are exhibited, respec- 
tively, in FIGS. 1 and 2. 

But how does a x scale as a function of t? When t is 
sufficiently small and the dispersion is large, one finds typi- 
cally that a x decreases as t~ x ^ 2 . This is consistent with the 
existence of only one population, but does not prove that 
multiple populations do not exist: A dispersion a x oc t~ p 
with p ~ 1/2 is also consistent with two distinct popula- 
tions which, however, are offset by such a small amount as to 
be completely indistinguishable. For larger sampling times, 
more variety is seen. In some cases, coinciding with ener- 
gies and potentials where trapping is at best infrequent, a x 
continues to scale as t~ ly/2 for larger t. However, when trap- 
ping is more important one finds invariably that, for a finite 
range of times, the dispersion decreases much more slowly 
than the i -1 ' 2 dependence expected for a single population. 
The dispersion of each separate population contributing to 
the total N[x(t)] may perhaps decrease roughly as i -1 ' 2 , but 
the composite a x for A r [x(^)] decreases much more slowly. 

This behaviour is illustrated in the upper two panels of 
FIG. 3 where, for In t between about 3 and 9, the total dis- 
persion decreases much more slowly. If, however, one probes 
somewhat longer intervals, the distinction between popu- 
lations becomes erased as a single chaotic orbit eventually 
transits from one 'type' of chaos to another. At sufficienly 
late times, there is only one population for which, as pre- 
dicted, a x decreases as t~ x ^ 2 . As illustrated in the two lower 
panels of FIG. 3, for potentials and energies where no evi- 
dence for multiple populations is observed, the dispersion is 
well fit by a £ -1//2 law throughout. 

The type of trapping that can arise and how it is man- 
ifested by the value of a short time Lyapunov exponent can 
be gauged from FIGS. 4 and 5, which were generated for 
a single orbit with energy E = 4 evolved in the potential 
(2) with a — b — 1. The initial condition was so chosen 
that, at the outset, the orbit is wildly chaotic. However, at 
a time somewhat larger than t = 15500 the orbit became 
trapped near a regular region where it remains stuck for an 
interval At ~ 1000 before again becoming untrapped, The 
left hand panels of FIG. 4 exhibit projections of the or- 
bit into the (x,y), (y,z), and (z,x) planes for the interval 
15700 < t < 16212, an interval At = 512 during which the 
orbit looks nearly regular. The right hand panels exhibit the 
same orbit for 15700 < t < 16724, an interval twice as long. 
It is clear that, during the second half of the longer interval, 
a significant qualitative change has occured: The orbit no 
longer manifests the reflection symmetries x — > —x, y — > — y, 
and z — > —z that were evident during the near-regular phase 
and, as viewed in the (z, x) projection, the orbit is no longer 
centrophobic. Indeed, if the orbit be integrated for a some- 
what longer time it becomes so wildly chaotic that the three 
different projections are almost indistinguishable visually. 

FIGURE 5, a time series of short time Lyapunov expo- 
nents for this orbit, was generated by computing x(^) f° r 
successive intervals St = 1 and, for each data point, perform- 
ing a box-car average over 256 adjacent times. At both early 
and late times x exhibits considerable variability. However, 
these changes are much more modest than what is observed 
during the orbit segment's near-regular phase, when x drops 
to much smaller values. When the exponent only assumes 



values larger than x ~ 0.10, the orbit appears visually to be 
wildly chaotic, but a drop to values smaller than x ~ 0.07 
signals a transition to a much more regular appearance. 

FIGURE 6 exhibits the analogue of FIG. 4 for a differ- 
ent initial condition, now evolved with energy E — 10 in the 
potential (2) with a = 1.6 and b = 0.4. 

Three conclusions seems inevitable: (1) As in two- 
dimensional potentials, it is possible for chaotic orbits to 
become trapped near regular regions for very long times, al- 
though they will eventually escape. (2) The presence of such 
trapping correlates with the possibility of more complicated 
phase space transport, indicative of the fact that a single 
connected phase space region divides into seemingly distinct 
populations. (3) The existence of de facto populations over 
some, but not all, finite intervals can be quantified numeri- 
cally in terms of the statistics of N[x(t)], the distribution of 
short time Lyapunov exponents. 



3 FIRST PASSAGE TIME EXPERIMENTS 

The experiments described here constitute a three-degree- 
of-freedom generalisation of two-degree-of-freedom first pas- 
sage time experiments in Pogorelov & Kandrup (1999). 
These involved identifying chaotic orbits which, in the ab- 
sence of any perturbations, remain stuck near regular islands 
for very long times, and determining how the introduction 
of weak perturbations, idealised as noise or periodic driv- 
ing, reduces the escape time. To assess the effects of noise 
of fixed amplitude and form, large numbers of noisy integra- 
tions of the same initial condition were performed, each with 
different pseudo-random seeds. The effects of periodic driv- 
ing with frequencies 'near' uj, involved instead large num- 
bers of integrations performed with slightly different fre- 
quency selected from some small interval Auj. The exper- 
iments with frequencies 1 < u) < 100 typically involved uni- 
formly sampling an interval Auj — 1. For experiments with 
100 < uj < 1000, the interval was Au> = 10; for those with 

0. 1 < uj < 1.0, Auj = 0.1. Each set of computations involved 
between 1000 and 5000 orbits. 

In the absence of any perturbations, determining pre- 
cisely the phase space regions corresponding to a 'sticky' 
chaotic orbit is possible, albeit very tedious (Contopoulos, 
private communication). However, the notion of a well de- 
fined boundary necessarily evaporates as soon as the orbits 
are perturbed: because of the perturbations, energy is no 
longer conserved, so that the effective phase space hyper- 
surface on which the orbit moves changes continually as the 
system evolves. For this reason, 'escape' was defined in a 
more practical fashion. Specifically, simple polynomial for- 
mulae were used to delineate approximately the configura- 
tion space region to which the orbit was originally confined, 
and the first escape time was identified as the first time the 
orbit leaves this region. That the escape criterion is reason- 
able was tested in two important ways: It was verified that, 
with or without perturbations, small changes in the bound- 
ary change the escape time only minimally; and that, at least 
in the absence of perturbations, 'escape' coincides with an 
abrupt increase in the short time Lyapunov exponent. 

In analysing the data, two simple diagnostics proved 
especially convenient: 

1. The time T(0.01) required for one percent of the orbits 
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in the computation to escape. As described below, escape 
does not begin immediately. Rather, there is typically an 
extended initial period, the duration ol which depends on 
the form and amplitude of the perturbation, during which 
there are no escapes. Escapes then turn on abruptly, at a 
time well characterised overall by T(0.01). (T(0.01) is less 
sensitive to statistical fluctuations than To , the time the first 
orbit escapes.) 

2. The initial escape rate A. In most, albeit not all, cases es- 
capes, once they begin, appear to sample a Poisson process, 
at least initially, so that N(t), the number of orbits not yet 
having escaped, decreases exponentially. 
A good representation for the data, used in the analysis, was 

N(t) = JVoexpj - A[t-T(0.01)]} for t > T(O.Ol). (6) 

The experiments with periodic driving entailed solving 
an evolution equation of the form 

d 2 r 

^2 = - W(r) + A sin ujt f (7) 

for variable A and 10 Those with noise involved solving 
Langevin equations (cf. Chandrasekhar 1943b, van Kampen 
1981) 

J = -W(r) - V v + F. (8) 

Here r\ — r](v) is the coefficient of dynamical friction and 
the stochastic force F is homogeneous Gaussian noise, with 
moments 

(F a (t))=0 (a,b = x,y,z) 
and 

(F a (ti)F b (t 2 )) = K(v,t 1 -t 2 )5 ab . (9) 
For delta-correlated white noise, 

jr(v,T) = 20ij(v)fo(T), (10) 

with O a characteristic 'temperature,' i.e., a typical energy 
for the internal degrees of freedom responsible for the noise. 
Coloured noise replaces 5d (t) by a function of finite width. 
Most of the integrations involved choosing ~ E. 

Internal irregularities should give rise to both friction 
and noise, so that it is not reasonable physically to consider 
one without the other. However, the relative importance of 
friction versus noise as a source of phase space transport was 
tested by comparing computations including both friction 
and noise with computations with the friction rjv turned 
off. The white noise simulations were performed using an 
algorithm developed by Greiner et al (1988) [see also Hon- 
erkamp (1994)]. The coloured noise simulations were imple- 
mented using a new algorithm developed by I. V. Pogorelov 
as part of his Ph. D. dissertation. The basic idea is outlined 
in Pogorelov & Kandrup (1999). 

3.1 The effects of periodic driving 

As noted already, escape is a two stage process. At very 
early times there are no escapes, and the only significant 
effect of the driving is to cause different integrations of the 
same initial condition with different frequencies to diverge 
exponentially inside the trapped region. Overall, for a fixed 



frequency interval Aui, the rms dispersion in the separation 
of different orbits scales as 

8r rma oc A exp(xt), (11) 

where A is the amplitude of the driving and \ is compara- 
ble to a typical value for the largest short time Lyapunov 
exponent. Only after the orbit ensemble has dispersed to 
probe significant portions of the trapped region and 5r rm s 
approaches some critical value do any escapes occur. These 
turn on quite abruptly, the interval during which (say) the 
first percent of the orbits escape being much shorter than the 
time before the first orbit escapes. That escapes begin when 
Sr rma approaches some critical value implies that the time 
at which escapes begin, as probed, e.g., by the one percent 
escape time T(0.01), scales logarithmically in amplitude. 

The earlier stages of the actual escape process, during 
which 90 percent or so of the orbits escape, are typically well 
represented by a Poisson process, with N(t), the fraction of 
the orbits that have not yet escaped, decreasing exponen- 
tially. This allows for the simple interpretation that, once 
the orbits have diverged to probe all the trapped regions, 
they can and will escape 'at random' as they find an ap- 
propriate exit channel. However, one finds typically that, 
at later times, escapes proceed more slowly, so that the de- 
crease in N(t) becomes subexponential. The reason for this 
is not completely clear, but it seems reasonable to suppose 
that some of the orbits that did not escape early on be- 
came trapped even more closely to the regular region, so 
that escape is more difficult than initially. (Because energy 
is no longer conserved, it is possible for an initial condi- 
tion which, in the absence of irregularities, corresponds to a 
chaotic orbit to move into a phase space region where, in the 
absence of the irregularities, it would be regular, at which 
point escape becomes exceedingly difficult.) For very small 
amplitudes A, the escape rate A is often relatively insensi- 
tive to A, which suggests that the driving does not facilitate 
escapes all that much. However, for larger amplitudes, say 
A > 10~ 6 -10~ one typically finds that A exhibits a roughly 
logarithmic dependence on A. Increasing the amplitude in- 
creases the escape rate. This indicates that the driving helps 
the orbits access escape channels by jiggling them about. 

Periodic driving has the strongest effect for frequencies 
lu ~ t~jj , which allow an efficient resonant coupling between 
the driving frequency and the natural frequencies of the un- 
perturbed orbit. However, periodic driving can still have a 
significant effect at much larger and much smaller frequen- 
cies, presuming by couplings through harmonics. Even fre- 
quencies as large as u — 1000 can significantly reduce the 
escape time. A detailed plot of T(0.01) or A as a function 
of lo exhibits considerable structure but, overall, the depen- 
dence on frequency is roughly logarithmic. 

FIGURE 7 exhibits representative plots of In N(t), gen- 
erated for different integrations of the same initial condition 
with E = 10 evolved in the potential (2) with a = 0.4 and 
b — 1.6. Each ensemble of 1000 orbits was generated by freez- 
ing the amplitude at a fixed value A and uniformly sampling 
the interval 5 < uj < 6 with 1000 different driving frequen- 
cies. The different curves correspond to different amplitudes. 
It is clear that, in each case, there is a finite initial interval 
during which there are no escapes, followed by an interval 
during which N(t) exhibits a roughly exponential decrease. 
The amplitude dependence of the escape process is illus- 
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trated in FIGS. 8 (a) and (b), which exhibit T(0.01) for the 
same initial condition for two different frequency intervals, 
2 < lo < 3 and 5 < to < 6. FIGS. 8 (c) and (d) exhibit 
the frequency dependence of T(0.01), as extracted from a 
collection of simulations with fixed amplitude A — 10~ 3 . 

3.2 The effects of white noise 

The effects of white noise are very similar to those of peri- 
odic driving, a fact that can be understood if one recognises 
that, in a real sense, white noise is an incoherent sum of os- 
cillations with all possible frequencies. Gaussian white noise, 
with random phases and a flat power spectrum, is equivalent 
mathematically to a superposition of periodic oscillations 
with all possible frequencies and all possible phases! 

Overall, the effects of noise are again manifested as a 
two-stage process: a phase during which different orbits in 
the ensemble - now different noisy realisations of the same 
initial condition evolved to sample the same random process 
- disperse within the confining boundary, followed by a pe- 
riod of escapes reasonably well approximated as a Poisson 
process. FIG. 9 illustrates the typical behaviour of lnjV(£) 
allowing for a constant coefficient of dynamical friction and 
additive noise, generated for the same initial condition used 
to create FIG. 7. The different curves have the same "tem- 
perature," O = E = 10.0, but different values of rj. As for 
the case of periodic driving, T(0.01) increases logarithmi- 
cally with decreasing amplitude - now measured by rj - at 
least for 77 not too small, say rj > 10~ 9 — 10~ 8 . The slope 
A associated with the exponential decay of N also tends to 
decrease with decreasing rj, which indicates that, even af- 
ter escapes have begun, noise accelerates the escape process 
by jiggling orbits and thus helping them to find an escape 
channel. FIGS. 10 (a) and (b) exhibit T(0.01) as a function 
of log rj for two different initial conditions. 

As for the case of periodic driving, the effects of the 
noise largely disappear when the amplitude becomes too 
small, rj < 10" 10 - 10" 9 or so. This can be understood at 
least in part as a numerical artifact. The white noise simu- 
lations were effected using a fixed St time step, fourth order 
Runge-Kutta integrator, which leads to "random" errors of 
order (St) 4 . Given a time step St = 10 -3 , the computa- 
tions might be expected to incorporate "numerical noise" 
of amplitude ~ 10 -10 , and additional irregularities of lower 
amplitude should have only minimal effects. 

Another significant conclusion is that both the quali- 
tative and quantitative effects of noise seem comparatively 
insensitive to the details. Allowing for at least some simple 
forms of multiplicative noise and/or allowing for a variable 
coefficient of dynamical friction or turning off the friction 
altogether has only a minimal effect. As an illustrative ex- 
ample, one can consider what happens if the quantity rj en- 
tering into both the friction and the autocorrelation function 
becomes a nontrivial function of velocity, assuming rj oc v p . 
The basic conclusion of such an investigation, illustrated 
by FIG. 11, is that these changes have almost no effect. 
Here the solid curve represents additive noise and a con- 
stant coefficient of dynamical friction, 770. The dot-dashed, 
triple-dot-dashed, and dotted curves represent, respectively, 
friction and noise with rj = rjov p /(v p ), with p = 1, 2, and 
3, and {v p ) a mean value computed for the unperturbed or- 
bit. (This normalisation ensures that the "average" noise 
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is the same for the additive and multipicative simulations.) 
The dashed curve corresponds to additive white noise but 
vanishing friction (i.e., allowing for a nonzero 770 in the au- 
tocorrelation function but assuming a vanishing coefficient 
of dynamical friction). The obvious point is that, early on, 
these curves are virtually identical and that, even at later 
times, the differences are comparatively minimal. 

In principle, one could perhaps artificially "tune" the 
form of the noise to enhance or suppress its effects, e.g., 
by making kicks especially large when an orbit approaches 
an escape channel. However, this seems unphysical. To the 
extent that the "noise" impacting real stars is largely un- 
corrected with the physics of the bulk potential, which de- 
termines the locations of these escape channels, one might 
anticipate that the details will be comparatively unimpor- 
tant. 

3.3 The effects of coloured noise 

The aim of the work described here was to determine how 
the aforementioned results regarding friction and white noise 
are altered if the noise becomes coloured. In other words, 
what happens if the autocorrelation function K(r) is not 
delta-correlated in time, so that the random impulses to 
which the orbits are subjected are of finite duration? 

Two comparatively simple examples were considered. 
The first corresponds to the so-called Ornstein-Uhlenbeck 
process (cf. van Kampen 1981), for which K(t) decays ex- 
ponentially: 

K(t) = aij6)exp(-a\T\). (12) 

The second alternative involved an exponential modulated 
by a power law: 

K( T ) = exp(-a|r|) ( 1 + a|r| + — r 2 V (13) 

8 V 3 / 

In each case, the normalisations were so chosen so that 

K{r)dT = 2Qrj. (14) 

- 00 

In other words, fixed noise amplitude means a fixed value 
2677 for the time integral of the autocorrelation function 
(which equals the diffusion constant D entering into a 
Fokker-Planck description [cf. Chandrasekhar 1943b, van 
Kampen 1981]). The autocorrelation times for these pro- 
cesses are, respectively, t c — l/a and t c = 2/a. The white 
noise calculations described in the preceding section can be 
understood as involving a singular limit a — *■ 00. As for the 
case of multiplicative noise, these two examples only probe 
the tip of an iceberg. However, an analysis of their effects 
does provide insight into the question of how a finite au- 
tocorrelation time can impact phase space transport in a 
complex phase space. 

As for white noise, the evolution of an ensemble of or- 
bits in the presence of coloured noise is a two stage process. 
After an initial interval without escapes, during which differ- 
ent members of the ensemble diverge exponentially, escapes 
turn on abruptly, with the first percent of the orbits escaping 
within a time T(0.01) much shorter than the time before the 
first escape. This is then followed by a phase during which 
N(t) decreases systematically. As for the case of white noise, 
this phase can usually be well fit overall by an exponential, 
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although a plot of N(t) for coloured noise tends to exhibit 
a bit more structure than do comparable plots for white 
noise or periodic driving. For fixed autocorrelation time tc, 
T(0.0f) exhibits a roughly logarithmic dependence on 77, at 
least for 10~ 8 < 77 < fCP 4 or so; and the slope A associ- 
ated with the exponential decrease in N(t) tends to scale 
logarithmically with 77. One also finds that the presence or 
absence of friction is largely immaterial. 

In all this, coloured noise behaves just like white noise. 
The real question is: how do things depend on i c ? When 
t c is very small and a is very large, the effects are nearly 
indistinguishable from white noise: both T(0.01) and A are 
essentially unchanged. Significant deviations only begin to 
arise when the autocorrelation time t c becomes comparable 
to the dynamical time tD for the unperturbed orbit. At this 
stage, the noise begins to have an appreciably weaker effect, 
its overall efficacy scaling logarithmically in a or t c . 

FIGURE 12 exhibits In N(t) for a fixed ensemble 
evolved with coloured noise satisfying eq. (13) for O = 10 
and r) = fCP 5 , but allowing for several different values of a. 
FIGURE 13 shows examples of how T(0.01) scales with 77 
for fixed nonzero a (a and c) and how T(0.01) scales with a 
for fixed 77 (b and d). 

The observed behaviour is easy to understand. Increas- 
ing t c from zero to a finite value is equivalent to replacing 
the flat white noise power spectrum by a band-limited power 
spectrum. Deviations from white noise begin to be impor- 
tant when this spectrum is so limited that there is compar- 
atively little power at frequencies comparable to the natural 
frequencies of the unperturbed orbit. In this sense, it would 
appear that, like modulational diffusion associated with pe- 
riodic driving, noise-induced extrinsic diffusion should be 
interpreted as a resonance phenomenon. The fact that the 
efficacy of the perturbation scales logarithmically in a is 
reminiscent of the fact that, for periodic driving, the effi- 
cacy scales logarithmically in driving frequency. 



4 CHAOTIC MIXING IN THE PRESENCE OF 
PERIODIC DRIVING AND NOISE 

The aim of the work summarised in this Section was to deter- 
mine how chaotic mixing is impacted by low amplitude per- 
turbations idealised as periodic driving or friction and noise. 
This entailed selecting localised ensembles of initial condi- 
tions and tracking their behaviour as they are integrated 
into the future, both with and without perturbations. The 
resulting orbital data were analysed as in Kandrup (1998) 
by computing both (i) coarse-grained representations of re- 
duced distribution functions f(Z a ,Zb) for different pairs of 
phase space coordinates and (ii) time-dependent moments 
(x i y 3 z k p l x p y n p z ) for i + j + k + I + m + n < 4. 

The principal conclusion is that low amplitude pertur- 
bations can impact the evolution in two potentially signifi- 
cant ways: 

Even in settings where "stickiness" is not important, so 
that orbits can spread out with comparatively few obstacles, 
time-dependent perturbations can play a role in damping 
oscillations and "fuzzing out" short wavelength structures. 
Because of the exponential sensitivity generic for chaotic sys- 
tems, a localised ensemble of initial conditions correspond- 
ing to "unconfined" chaotic orbits will, in the absence of 



any perturbations, diverge exponentially. However, this ten- 
dency of nearby orbits to diverge implies that the dispersing 
ensemble will eventually converge towards an invariant, or 
near- invariant, distribution, i.e., a (near-)equilibrium (Kan- 
drup 1998; see also Merritt & Valluri 1996). If, e.g., the or- 
bital data be binned at succesive instants to generate a grid- 
ded representation of some f(Z a , Zb,t), one finds that, with 
respect to a discrete L 1 or L 2 norm, this time-dependent 
f(t) will typically evolve exponentially towards some nearly 
time-independent f„i v (Z a , Zb), i.e., 

{Y, J2 \K Za > Zb > *) - fniv(Za,Z b )\^ ^ ~ e~ At , (15) 
a b 

for p — 1 or 2. Similarly, moments like (x) or {xp y ), which 
vanish for an invariant distribution, evolve towards zero ex- 
ponentially. 

However, this evolution is not always uniform in time. 
In many cases, the systematic evolution is accompanied by 
coherent oscillations that only damp on comparatively long 
time scales. For example, a plot of the dispersion in some 
phase space variable can exhibit significant oscillations that 
persist for tens of dynamical time to even if the "average" 
dispersion has settled very nearly to its equilibrium value 
within a time t < lOtzj. Allowing for noise and/or periodic 
driving can significantly decrease the amplitude of these os- 
cillations. 

This behaviour is illustrated in FIG. 14 (a) and (c) 
which track a x for the same ensemble of initial conditions 
evolved in the potential (b) with a = b = 1 both with and 
without periodic driving. In each cases, the dispersion con- 
verges towards a value a x w 1.2, but for the perturbed orbits 
the convergence is clearly more efficient. FIGS. 15 (a) and 
(c) track <j x for the same ensemble, now evolved in the pres- 
ence of additive white noise with = E and, respectively, 
77 = 10 -4 and 77 = 10 -6 . It is clear that weak noise can 
damp irregularities even more efficiently than can periodic 
driving. 

The second important point is that, because of the es- 
cape phenomenon described in the preceding section, weak 
perturbations can accelerate the diffusion of sticky orbits, 
allowing them to probe phase space regions which otherwise 
are only accessed at comparatively late times. For exam- 
ple, as noted in Kandrup (1998), ensembles of chaotic orbits 
with small initial z and p z tend, in the absence of perturba- 
tions, to remain relatively close to the z-axis, so that they 
only approach an invariant distribution on a relatively long 
time scale - t ~ 50 rather than t ~ 5. However, low ampli- 
tude perturbations can reduce this time scale significantly 
by allowing these orbits to move away from the z-axis, thus 
enabling them to probe all (or at least a signficantly larger 
portion of) the connected phase space region. 

This is illustrated in FIGS. 14 (b) and (d), which track 
cr 2 for the same ensemble used to generate FICS. 14 (a) and 
(c). Here once again the top panel represents an unperturbed 
orbital integration whereas the lower panel is generated from 
a simulation subjected to a finite amplitude periodic driving. 
In interpreting the plots of a x and a z , it should be noted 
that the initial ensemble was chosen so that x — 0.0, z = 
0.2, and p z — 0.0. The initial y and p y were generated by 
uniformly sampling, respectively, the intervals [0.8, 1.0] and 
[2.4, 2.6] and p x > was computed as a function of the 
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remaining phase space variables to ensure that E = 10. The 
unperturbed dispersion a z in FIG. 14 (b) is typical of what 
is seen when evolving a collection of relatively sticky orbits. 
FIG. 14 (d) more closely resembles what one would have 
expected to see for an unperturbed evolution of orbits that 
are not initially trapped in a special phase space region. 
FIGS. 15 (b) and (d) exhibit analogous plots of a z generated 
from the noisy simulations used to generate FIGS. 15 (a) and 
(c). The trend is again the same, although the effect is not 
quite as pronounced 

What this could imply for the visual appearance of a 
galaxy is illustrated in FIGS. 16 and 17, which exhibits grey 
scale plots generated from the same orbit integrations. The 
first panel in FIG. 16 exhibits the distribution f(z,y) for 
an unperturbed integration, generated from data recorded 
at intervals St = 1 for 64 < t < 128. The obvious point is 
that, even though the potential exhibits cubic symmetry, 
the orbits are more localised in the z-direction than in the 
y- (or as—) direction. The second panel exhibits f(z, y) much 
later on, namely for 896 < t < 1024, by which time the dis- 
tribution has become much more symmetric. The lower pan- 
els, again generated for 64 < t < 128, exhibit f(z, y) for the 
same initial ensemble now evolved in the presence of periodic 
driving with two different frequencies. FIGS. 17 (a) and (b) 
exhibit analogous plots generated from noisy integrations, 
It is evident in each case that, because of the perturbations, 
the distribution is more nearly symmetric than was the case 
for an integration with no perturbations. 



5 DISCUSSION 

To assess the importance of all this for galactic dynamics, 
one must determine the extent to which orbit trapping and 
slow phase space diffusion, seemingly generic phenomena for 
nonintegrable potentials that admit a coexistence of both 
regular and chaotic orbits, are present in the cuspy, traxial 
potentials which appear to characterise many galaxies. 

As a concrete example, consider the maximally triaxial 
7 = 1 Dehnen potential (Merritt and Fridman 1996), which 
is generated self-consistently from the mass density 



(3-t)M 
p{m) — —. — i- — m 



with 



4mabc 



2 2 2 

x y z 
a 2 b 2 c 2 ' 



'(l + s 



-(4- 7 ) 



(16) 



(17) 



for 7 = 1, c/a = 1/2, and (a 2 - b 2 )/{a 2 - c 2 ) = 1/2. For 
each of twenty different energies, some 300 — 400 well sep- 
arated initial conditions were identified, all corresponding 
to chaotic orbits. [These were the chaotic initial conditions 
used to generate Siopis's (1998) library of orbits for the con- 
struction of equilibria using Schwarzschild's (1979) method.] 
Each initial conditions was integrated for a time 2048£d, 
with x{t) recorded at intervals St = 1. As for the compu- 
tations described in Section 2, the resulting data were then 
partitioned into segments of various length t and the disper- 
sion a x was studied as a function of t. 

The principal conclusion from this analysis is that, at 
least for this potential, unperturbed phase space transport 
always proceeds very slowly. For relatively high energies, the 
observed distributions of short time Lyapunov exponents are 



similar to what is found for the potential (2) for choices of 
the parameters a and b and the energy E for which "stick- 
iness" is a frequent occurence. Moreover, as for that poten- 
tial, the time scale on which details wash out and the disper- 
sion begins to decay as t^ 1 ^ 2 is of order lOOOiij. However, 
for the lowest energies, where chaotic orbits pass very close 
to the central cusp, phase space transport appears to pro- 
ceed even more slowly! Distributions of short time Lyapunov 
exponents tend to exhibit considerable structure, even for 
comparatively long times. Indeed, even after 2048to there 
is no suggestion that the distinctions between different pop- 
ulations have been erased. This is illustrated in FIG. 18, 
which exhibits a x (t) for the two highest and lowest energies 
used by Merritt & Fridman (1996) in their construction of 
Schwarzschild equilibria. 

The obvious inference is that, at least for this triaxial 
Dehnen potential, unperturbed chaotic orbits often diffuse 
through phase space only very slowly, especially at low en- 
ergies. The two obvious questions then are: (i) is this result 
generic for cuspy triaxial potentials, or just an accident for 
this specific potential, and (ii) to what extent can low ampli- 
tude friction and noise accelerate even this very slow phase 
space transport? Both these questions are currently being 
addressed (Siopis & Kandrup 1999). 
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Figure 1. N[\(t)], the distribution of short time Lyapunov expo- 
nents for a collection of chaotic orbits with energy E = 6 evolved 
in the potential (2) with a = 3 and b = — 1, computed as a func- 
tion of the sampling time t. (a) t = 1. (b) t = 4. (c) t = 16. (d) 
t = 64. (e) t = 256. (f) t = 1024. (g) t = 4096 (h) t = 16384. 
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Figure 3. (a) The dispersion <r x (t), plotted as a function of 
the sampling time t, for a collection of chaotic orbits with energy 
E = 4 evolved in the potential (2) with a = 6 = 1. The dashed 
line corresponds to a t -1 / 2 dependence, (b) The same for chaotic 
orbits with E = 6 and a = b = 1. (c) The same for chaotic orbits 
with E = 6 and a = 3 and fe = —1. (d) The same for chaotic 
orbits with E = 10 and a = 3 and fe = —1. 



Figure 2. TV [x(*)] > the distribution of short time Lyapunov expo- 
nents for a collection of chaotic orbits with energy E = 6 evolved 
in the potential (2) with a = b = 1, computed as a function of 
the sampling time t. (a) t = 1. (b) t = 4. (c) t = 16. (d) t = 64. 
(c) t = 256. (f) t = 1024. (g) t = 4096 (h) t = 16384. 



[t] 



[t] 




Figure 4. Visual representations of the orbit used to compute 
Fig. 4. The left hand panels exhibit data for times 15700 < t < 
16212. The right hand panels exhibit data for 15700 < t < 16724, 
an interval twice as long. 



Figure 6. The analogue of Fig. 5 for an initial condition with 
E = 10, in this case integrated in the potential (3) with a = 1.6 
and b = 0.4. In this case, the integration times for the left and 
right hand panels are At = 400 and At = 800. 
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Figure 5. A coarse-grained estimate of the largest short Lya- 
punov exponent \, generated for a chaotic orbit in the potential 
(3) with E = 4 and a = b = 1. The coarse-graining was effected 
by computing x(t) f° r successive intervals St = 1, and performing 
a box-car average over 256 adjacent intervals. 




Figure 7. N(t), the fraction of the orbits from a 1000 orbit 
ensemble not yet having escaped at time t, computed for an initial 
condition with E = 10.0 evolved in the dihedral potential with 
a = 0.4 and b = 1.6 in the presence of periodic driving with 
5 < lu < 6 and variable amplitude A = 10~ 2 (solid line), A = 
10~ 3 (dashed), A = 10~ 4 (dot-dashed), A = 10~ 5 (triple-dot- 
dashed), and A = 10~ 7 (dotted). 




Figure 8. (a) T(0.01), the first escape time for 1% of an ensemble 
of 2000 integrations of the initial condition in Fig. 6, driven with 
frequencies 2 < u> < 3 and variable amplitude, (b) The same for 
5 < u> < 6. (c) T(0.01) for the same initial condition, now plotted 
as a function of ui for fixed amplitude A = 10~ 3 . (d) T(0.01) 
plotted as a function of log ui. 
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Figure 10. (a) T(0.01), the first escape time for 1% of an ensem- 
ble of 2000 integrations of the initial condition in Fig. 6, perturbed 
by additive white noise with © = 10.0 and variable rj. (b) The 
same for a different initial condition. 
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Figure 9. N(t), the fraction of the orbits from a 1000 orbit 
ensemble not yet having escaped at time t, computed for an initial 
condition with E = 10.0 evolved in the dihedral potential with 
a = 0.4 and b = 1.6 in the presence of additive white noise with 
© = 10.0 and variable rj = 10~ 3 (solid line), r\ = 10~ 4 (broad 
dashes), rj = 10 -4,5 (short dashes), rj = 10~ 5 (dot-dashed), rj = 
10~ 5 ' 5 (triple-dot-dashed), and rj = 10" 8 (dotted). 



Figure 11. N(t), the fraction of the orbits from a 1000 orbit 
ensemble not yet having escaped at time t, computed for an initial 
condition with E = 10.0 evolved in the dihedral potential with 
a = 0.4 and b = 1.6 in the presence of different white noises with 
= 10.0 and rjo = 10~ 4 . The solid curve represents additive 
white noise and friction with constant rj. The dot-dashed, triplc- 
dot-dashed, and dotted curves represent, respectively, friction and 
noise with rj oc v, v 2 and v s . The dashed curve represents additive 
white noise but vanishing friction. 
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Figure 12. N(t), the fraction of the orbits from a 4800 orbit 
ensemble not yet having escaped at time t, computed for the 
initial condition exhibited in Fig. 9, allowing for friction and col- 
ored noise given by eq. (5.2) with © = 10, r\ = 10~ 5 , and either 
white noise (triplc-dot-dashed) or variable a = 2.0 (dot-dashed), 
a = 0.2 (broad dashes), and a = 0.02 (solid). 
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Figure 13. (a) T(0.01), the first escape time for 1% of an ensem- 
ble of 4800 integrations, computed for the initial condition used 
to generate Fig. 9, plotted as a function of log?; for fixed a = 2.0 
for the stochastic processes denned by (5.1) (diamonds) and (5.2) 
(triangles), allowing for both friction and noise, (b) T(0.01) for 
the same initial condition, plotted as a function of log a for fixed 
r\ = 10~ 5 , for the stochastic processes (5.1) (diamonds) and (5.2) 
(triangles), again allowing for both friction and noise. The dashed 
line represents the asymptotic value for white noise (a — > oo). (c) 
The same as (a), albeit for a different initial condition and with 
a = 0.2. (d) The same as (b), albeit for the initial condition in (c) 
and with rj = 10~ 7 . as (a) and (b) for another initial condition. 
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Figure 14. (a) The dispersion a x generated for a localised en- 
semble of 6400 initial conditions with E = 10 evolved without 
perturbations in the potential (2) with a = 1 = 1. (b) The dis- 
persion a z for the same ensemble, (c) a x for the same ensemble, 
now allowing for periodic driving with £1 = 1 and A = 10 _l s . 
(d) a x for the same ensemble, allowing for periodic driving with 
$1 = 1 and A = lO" 15 . 




Figure 15. (a) The dispersion cr x generated for the same ensem- 
ble evolved in the presence of friction and additive white noise 
with = 10 and r\ = 10~ 4 . (b) The dispersion a z for the same 
computation, (c) The same as (a) but with r\ = 10~ 6 . (d) The 
same as (b) but with r\ = 10~ 6 . 
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Figure 16. (a) Gray scale plot of the distribution f(x, y) gen- 
erated for a localised ensemble of 6400 initial conditions with 
E = 10 evolved without perturbations in the potential (2) with 
a = b = 1, with data recorded at intervals St = 1 for 64 < t < 128. 
(b) The invariant distribution associated with a much longer time 
integration, (c) The same as (a) but allowing for periodic driv- 
ing with n = 1 and A = 10~ 1,5 . (d) The same as (c) but with 
A= 1Q-' 2 - 5 . 
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Figure 17. (a) Gray scale plot of the distribution f(x, y) gen- 
erated for the same initial ensembles as in the preceding figure, 
now allowing for additive white noise with = 10 and r; = 10~ 4 . 
(b) The same as (a) but with r\ = 10~ 6 . 
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Figure 18. The analogue of Fig. 3 for chaotic orbits in the max- 
imally traxial 7 = 1 Dchncn potential. The upper panels are for 
the two highest energies used by Merritt & Fridman (1996). The 
lower panels are for the two lowest energies. 



